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Abstract 



> 

o 

o 

^^ I We consider the problem of separation of variables in the Kramers 

Qs. ' equation admitting a non-trivial symmetry group. Provided the ex- 

0\ • ternal potential V{x) is at most quadratic, a complete solution of the 

r~| [ problem of separation of variables is obtained. Furthermore, we con- 

C^' struct solutions of the Kramers equation with separated variables in 

r^ ■ explicit form. 

B 

t^- ' 1 Introduction 

> 

/\ • Many phenomena in physics and, especially, in chemical physics may be mod- 

j^ ■ elled as the Brownian motion of particles in an external potential V{x), the 

appropriate transport equation being the (l-|-2)-dimensional Fokker-Plank 
equation of special form 

ut = vuyy - yu^ + {vy + V'{x))uy + z/M, (1) 

where u = u(t,x,y) is a sufficiently smooth real-valued function and j/ is a 
real parameter. 

The first relevant result on studying partial differential equation (PDE) 
(|I]) has been obtained by Kramers |]1|. He found a solution of the escape 
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problem of a classical particle subjected to Gaussian white noise out of a deep 
potential well. This is why the equation in question is called the Kramers 
equation (KE) (see, for more details 0]-[0). 

As KE is a PDE with variable coefficients, we cannot apply the Fourier 
transform in order to solve it. In fact the only way to obtain exact solutions of 
KE are either to utilize its Lie symmetry or to apply the method of separation 
of variables. The first possibility has been exploited recently in |^ §, where 
symmetry classification of the class of PDEs (||) has been carried out. The 
principal result of these papers is that KE has a symmetry group that is 
wider than a trivial one-parameter group of time translations if and only if 
V"{x) = 0. 

The principal aim of the present paper is to apply the direct approach 
to variable separation in PDEs suggested in 0^0 to solve KE. As is well- 
known, separability of PDE is intimately connected to its symmetry within 
the class of second-order differential operators [|T0|. This is why, we will 
concentrate on the case V{x) = kx, k = constant, namely, we consider KE 
having non-trivial Lie symmetry 

ut = vuyy — yux + iyy + kx)uy + vu. (2) 

In a classical setting the method of separation of variables (say, in the 
Cartesian coordinate system) is based on a special representation of a solution 
to be found in factorized form: 

u{t,x,y) = (poit)(piix)(p2iy), 

where (fi,i = 0,1,2 are solutions of some ordinary differential equations 
(ODEs). However, one can try to separate variables in the equation under 
study in another coordinate system, for example, in polar coordinates and 
look for a solution of the form 



n(t, X, y) = ipoit)(fi Ux'^ + yA ip2 i 



y 

arctan — 

X 



So, if we are given some coordinate system, then it is clear how to get exact 
solutions with separated variables. However, the classical approach gives no 
general routine for finding all possible coordinate systems providing sepa- 
rability of the equation under study. Our approach to the problem of the 
separation of variables in evolution-type equations (to be specific, we take 



the case of an equation having three independent variables t, x, y) is based 
on the following observations: 

• All solutions with separated variables known to us can be represented 
in the form 



u{t,x,y) = Q{t,x,y)ipo{t)(fi{uJi{t,x,y))ip2{uj2it,x,y)), 



(3) 



where Q, ui, UO2 are sufficiently smooth functions and ipi,i = 0, 1, 2 sat- 
isfy some ODEs. 

• The functions (f i,i = 0,1,2 depend on two arbitrary parameters Ai, A2 
called spectral parameters or separation constants. Furthermore, the 
functions Q,uJi,U2 are independent of Ai, A2. 

By properly postulating these features we have formulated an efficient 
approach to the problem of variable separation in linear PDEs [^ . Applying 
it to the KE (|]) we look for its particular solutions of the form (^, where 
functions Q, <jJ\,ijJ2 are chosen in such a way that inserting (§) into KE yields 
three ODEs for functions y9o(^), V'i(^i), V'2(^2) 



f/o(i,V'o,0o;Ai,A2) = 0, 
\Ji{u)i, ipi, (pi, ipi, Ai, A2) = 0, 



1,2. 



(4) 



Here Uq, Ui, U2 are some smooth functions of the indicated variables, Ai, A2 
are real parameters and, what is more. 



rank 
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2. 



(5) 



Note that the functions Q,uji,uj2 are independent of Ai, A2. 

Provided the above requirements are met, we say that KE is separable in 
the coordinate system t,uji{t,x,y),uj2{t,x,y). 

Due to the fact that the equation under study is hnear, the reduced 
equations prove to be linear as well. Furthermore, we have to consider two 
distinct cases. 



Case 1. The system of equations (^ has the form 

00 = Ao{t; Xi, X2)y:>o, 

ipi = Ai{uJi; \i, \2)'fi, (6) 

(p2 = A2{UJ2] Xl, X2)^2 + A3{UJ2] Xl, X2)<^2- 

Case 2. The system of equations @ has the form 

00 = ^o(i;Ai,A2)v5o, 

01 = Ai(u;i; Ai,A2)v5i, (7) 

In these formulae Aq, . . . , A3 are some smooth real- valued functions of the 
indicated variables. 

Consequently, there are two different possibilities to separate variables in 
KE, either to reduce it to two first-order and one second-order ODEs or to 
three first-order ODEs. It is impossible to reduce KE to two or three second- 
order ODEs because it contains a second-order derivative with respect to one 
variable only. 

Provided the system of reduced ODEs has the form (j^), separation of 
variables in (^ is performed in the following way: 

1. We insert the Ansatz (|) into KE and express the derivatives 0o, 0i, 
01, 02 in terms of functions ipo, ipi, ip2, 02 using equations ® and their 
differential consequences (where necessary). 

2. The equality obtained is split by ipo, y^i, ^2-, 02, Ai, A2 which are regarded 
as independent variables. This yields an over-determined system of 
nonlinear PDEs for unknown functions Q^lli\^llI2- 

3. After solving the above system we get an exhaustive description of 
coordinate systems providing separability of KE. 

Clearly, if we adopt a more general definition of the separation of vari- 
ables, then additional coordinate systems providing separability of KE may 
appear. However, all solutions with separated variables of the Schrodinger 
and heat conductivity equations known to us can be obtained within the 
described approach. 



The case when the system of reduced ODEs is of the form (|^ is handled 
in a similar way. 

Next, we introduce an equivalence relation on the set of all coordinate 
systems providing separability of KE. We say that two coordinate systems 
t, ui, uji and t', cu^, uj^ are equivalent if the corresponding solutions with 
separated variables are transformed one into another by 

• the group transformations from the Lie transformation group admitted 
byKE, 

• the transformations of the form 

t^t' = fo{t), uJi-^uj- = fi{uji), (8) 

Q^Q' = Qho{t)h{u;,)h2{u;2), (9) 

where fo, fi,ho, hi are some smooth functions. 

It can be proved that formulae (§), (P) define the most general transfor- 
mation preserving the class of Ansatze (^. The equivalence relation split 
the set of all possible coordinate systems into equivalence classes. In a se- 
quel, when presenting the lists of coordinate systems enabling us to separate 
variables in KE we will give only one representative for each equivalence 
class. 

2 Principal results 

In this section we give a complete account of our results on the separation 
of variables in KE obtained within the framework of the approach described 
in Introduction. We write down explicit forms of the functions Q{t,x,y), 
uji{t,x,y), U2it,x,y) and the corresponding reduced ODEs for functions 

ipo{t), V2i(c^i), ^2iuJ2). 

Theorem 1 Equation ^ admits the separation of variables into two first- 
order and one second-order ODEs if and only ifk takes one of the three values 
0, 3z/^/16, — 3z/^/4. Furthermore, equations separates into three first-order 
ODEs with an arbitrary k. 

Theorem 1 gives a general description of separable KEs. The solution of 
the problem of separation of variables in corresponding KEs is provided by 
Theorems 2-6 later. 



Theorem 2 The set of inequivalent coordinate systems providing separabil- 
ity of KE with k = i/^/4 is exhausted by the following ones: 

^ Jiy_=Lli^ 2 = 1,2, 
/2/1-/1/2' 



(10) 



+ |-&^^^^^--U^--ln|/2/i-/i/2| + -t|; 

\^^f2fl-flf2 4; 2 '^^'^ '^''' 2/' 

(fj^i + h^V . . . . 

^o = ^\ -y-. Yj- ^0, ^1 = Ai^i, v?2 = A2^2, 

V /2/I - Jl J2 / 

where 

fi = t i Ai sinh —t + Ao cosh —t] + A-i sinh —t + Aa cosh — t, 
Jl \^i 2 2/ 2 2 

/2 = t ( -Bi sinh —t + i?2 cosh — t 1 + -B3 sinh — t + i?4 cosh — t, 

and Ai,...,B4 are arbitrary real constants satisfying the condition 2Ci2 
—^{Ci3 — 6*24)= 0. Hereafter we use the notations 

Qj = BiAj - AiBj, i,j = l,...,4. 

Theorem 3 The set of inequivalent coordinate systems providing separabil- 
ity of KE with k > z/^/4 is exhausted by those given in ( jT^ j with 

/i = sin bt{Ai sinh at + A2 cosh at) + cos bt{As sinh at + A4 cosh at), 

/2 = sin bt{Bi sinh at + B2 cosh at) + cos bt{B3 sinh at + -B4 cosh at). 



where a = ^,b = Jk — ^ and Ai, . . . , S4 are constants fulfilling the condition 
(C12 + 6*34)6 + (Ci3 — C24)a = 0. The explicit form of the function Q and 
the reduced ODEs are also obtained from the formulae ( |7Z[ ) with /i, /2 given 
previously. 



Theorem 4 The set of inequivalent coordinate systems providing separabil- 
ity of KE with k < z/^/4 and k ^ 0, 3i/^/16, — 3z/^/4 is exhausted by those 
given in ( |T^ with 

/i = sinh bt{Ai sinh at + A2 cosh at) + cosh btlA^ sinh at + A4 cosh at), 

/2 = sinh bt{Bi sinh at + B2 cosh at) + cosh bt{B3 sinh at + i?4 cosh at), 



where a = f,&= y^ — k and Ai, . . . , i?4 are constants fulfilling the condition 
(C12 — Cu)b + (C13 — C24)a = 0. The explicit form of the function Q and 
reduced ODEs are also obtained from the formulae ( |7^ j with /i, /2 given 
before. 



Theorem 5 The set of inequivalent coordinate systems providing separabil- 
ity of KE with k = is exhausted by 

1) those given in ([7^ with 

/i = Ai sinh vt + A2 cosh vt + A^t + A^, 

/2 = Bi sinh vt + i?2 cosh vt + i?3t + i?4, 

where Ai, . . . , S4 are constants fulfilling the equation vCi2 — C34 = 0. 
The explicit form of the function Q and reduced ODEs are also obtained 
from the formulae (fJD) with /i, /2 given previously; 



2) the following coordinate system: 



c^i = x, uj2 = y, Q = exp I - — 



ipO = l^\lLpO, (^1=1^X2^1, V32 = f — + A22/ + Ai - - j (y92. 



Theorem 6 The set of inequivalent coordinate systems providing separabil- 
ity of KE with k = 3z/^/16 or k = —Su^/A is exhausted by 

1) those given in Theorem 4 under k = 3z/^/16 or k = — 3z/^/4; 



2) the following coordinate systems: 
uJi = R^x, UJ2 = Ry + 3i?a;, 



3 R 15RR 



AuR 4i/i?2 
-^]x' + ^t + 2\nR\; 

ipo = jyXiR^^o, 01 = i^A2V5i, ^2 = {X2UJ2 + Xi)ip2, 



where 



m 



1 

cosh at ' 

1 
sinh at ' 

exp{±at} 



with a 



4' 


under k = 


3u' 
16' 


2' 


under k = 


3z/2 

4 



3 Proof of Theorems 1—6 

In order to prove the assertions of the previous section one should apply to 
equation (^ the algorithm of variable separation described in Introduction. 
We give a detailed proof for the case when system of reduced ODEs is of 
the form (j^). Inserting Ansatz (^ into KE (§) and expressing the derivatives 
00, 01) 01) 02 in terms of functions (/jq, V^i, V^2, 02 with the use of equations 
(P) and their differential consequences yield a system of two nonlinear PDEs 

QuJ2t + yQuJ2x = ^{yQ^2y + 2QyUJ2y + 2QAiUJiyUJ2y 

+QA2Ujly + QuJ2yy) + kxQu2y, (11) 

Qt + QAq + QAiUu + yQx + yQAic^i^. = v{Q + yQ^ 

+QAicuiyj/ + QA^uoly) + A;x(Qy + QAi^iy). (12) 

This system is to be split with respect to variables Ai, A2 (we remind that 
the functions uJi,U2 are independent of Ai, A2). To this end we differentiate 



11) with respect to Aj and get the relation 



{2Aix^iJly + A2X^U]2y)^2y = 0, 



1,2 



Due to the fact that uj2y does not vanish identically (otherwise it follows from 
( p!ID that UJ2 =constant), the equation 

2Aix,UJiy + A2X^UJ2y = 0, 2=1,2. (13) 

holds. 

Let us show first that we can, without loss of generality, put ooiy = 0. 
Suppose the inverse, namely that the inequality uiy ^ holds true. It 
follows from the second equation of system (^ that A3;^\^ + ^ia2 t^ 0- Let 
the function A\\^ be non-vanishing, then by the influence of (|13]) A2\^ ^ 0. 
Denoting 

^lAi = 5'(^i, Ai, A2), -2A2A1 = /(^2, Ai, A2) 



(14) 



we rewrite (^3]) as follows 

^ly /('^2,Ai,A2 



UJ2y 5f(t^l, Al, A2) 

Differentiating ([T^) with respect to Ai yields 

/ai _ 5'Ai 
/ " ^ ■ 

Hence we conclude that there is a function /c = /i;(Ai, A2) such that 

^ = ^ = A;(Ai,A2). 

Integrating the above equations we get 

/ = fci(Ai, A2)/i(c^2, A2), g = ki{Xi, A2)fi'i(c^i, Ai), 
so that (0) reduces to the relation 

^ly _ /l(^2,A2) 
^2y gi{uJl,X2) 

In a similar way we establish that the last relation is equivalent to the fol- 
lowing one 

^ly _ f2{uJ2) 

9 



hence 

Taking into account the equivalence relation (H) we can put g2 = 1 and 
/2 = 1 in the above equality thus getting uiy = U2y. Integrating this PDE 
yields 

CiJi = (^2 + h(t, x) 

with an arbitrary smooth function h. In view of this equation relation (0) 
takes the form 

2^A, +^A, =0, i = l,2. 

Hence we conclude that there exists a function A(Ai, A2) such that 

Ai = A(Ai,A2) + ii(a;i), A2 = -2A(Ai,A2) + i2(^2). (15) 

Within the equivalence transformation (^ with properly chosen functions 
hi,h2 we can put Ai{ui) = 0,A2{uj2) = 0. Furthermore, defining new sepa- 
ration constants as 

a; = a(Ai,A2), a'2 = A2 

and omitting the primes we represent (|15]) in the form 

Ai = Ai, A2 = -2Ai. 

Consequently, system (j^) takes the form 

V9o(t) = Ao{t)^o{t), 
0i(u;i) = Xiipi{uji), (16) 

V32(u;2) = -2AivJ'2(t^2) + ^3(^2, Ai, A2)v22(t^2)- 

Making the change of variables ip2 = 0exp{— AiCi;2} reduces the third 
equation of system (|16|) to 

= (A? + ^3)0. 

Let = 4>{uj2, Ai, A2) be a solution of this equation. Then the corresponding 
solution with separated variables becomes 



u = Q{t, X, y)(j){uj2, Ai, A2) exp <^ / Ao{t)dt + Xi{uJi - u;2) [> . 

10 



The structure of so obtained solution with separated variables is such that 
dependence of ui on y is inessential. Indeed, the function ui enters into the 
solution only as a combination ui — uj2 and the latter is equal to h{t,x). 
Consequently, we have proved that without loss of generality we may choose 

iOly = 0. 



Given the condition tuiy = 0, equation (|T3D reduces to the relations A2\^ = 
0, i = 1,2, hence we get A2 = A2{u2). Choosing appropriately the function 
/12 in d^) we can put A2 = 0. Next, differentiating (|^) with respect to Xi we 
arrive at the equations 

AoA, + AixX^it + y^ix) = i^A3x,uj2y, i = 1,2. (17) 

Differentiating twice the above equations with respect to y yields 

A3i^2UJ2Xi^2y + ^^3uJ2\i^2y^2yy + 2A3Xi{^2yy + ^2y^2yyy) = 0, (18) 

where i = 1,2. 

Note that due to ([17|) the inequality A^x. ^ holds. Dividing (|T^) into 
A^x^ and differentiating the equality obtained by Aj, j = 1, 2 we get 






3A, / x^ \ ^3K ) X, 

Case 1. At least one of the four expressions 

/ Az^2\i \ 



^2, + 5H^^ ^2yy = ^, I, J = 1,2. (19) 



A 



3A, / A 



does not vanish. Then it is easy to become convinced that the relation 

UJ2yy 



^2y 



fM 



holds true. Integration of this relation yields 

^2y = 9i {t, x) exp I / f{uJ2)duJ2 
where gi{t, x) is an arbitrary smooth function. 

11 



Next, by using the equivalence relation (|) we reduce the equation ob- 
tained to the form 

hence 

uJ2 = y9iit,x)+ g2{t,x), (20) 

g2{t,x) being an arbitrary smooth function. 

In view of this result, ([TsD takes the form As^^^jAi =0,2 = 1,2, hence 

As = Ai(Ai, X2)uj2 + A2(Ai, A2) + F{iU2), (21) 

where Ai,A2,-F are arbitrary smooth functions of the indicated variables. 
Furthermore, it is not difficult to prove that Ai, A2 are functionally indepen- 
dent (since otherwise the condition (^ would be broken) and, consequently, 
after redefining Ai, A2 we can represent (|2TD in the form 



^3 = AiCU2 + A2 + F(cU2). (22) 

Case 2. Suppose that now 



^3A. 



I =0, z,j = l,2. 

A, 



Integrating the above system of PDEs gives the following form of A^x. 

Asa, = Bi{uj2)Li{Xi, A2), i = 1, 2, (23) 

where Bi, Li are arbitrary smooth functions and, what is more, Bf + B2 7^ 0. 
As a compatibility condition of system (^) we get 

BlLix.^ = _B2-^2Ai- 

Subcase 2.1. L1X2 7^ 0, L2A1 7^ 0. Given this restrictions the compatibility 
condition is transformed to 

Bi{uJ2) L2X1 . ^„.s 

Integrating system (p3D with account of (p^ yields 

A3 = A(Ai,A2)Fi(cU2)+F2(cU2), 

12 



where A, F are arbitrary smooth functions of the indicated variables. After 
redefining separation parameters Ai, A2 we represent the relation as follows 

A3 = AiFi(^2)+F2(^2). (25) 



Subcase 2.2 Lw^ = 0, L2A1 = 0. Integrating system ( psf ) and redefining 
the separation parameters Ai, A2 yield 

A'i = XiSi{uj2) + A2S'2(u;2) + 5*0(^2), (26) 

where 5*1, 5*2, Sq are arbitrary smooth functions. An analysis of formulae (]2^), 



(pSD and (p6|) shows that the first two of them are particular cases of formula 
(P^D. Thus, the most general form of the function A3 is given by (PB|). 

Inserting ( ]26|) into (p!7D and differentiating the equality obtained with 
respect to x and A^ gives Aix-x^ = 0,i,j = 1,2. Hence, we get for Ai 

A, = AiLi(cji) + A2L2(c^i) + Lo(c^i), (27) 

where Li,L2,Lo are arbitrary smooth functions. 

Next, inserting (pB]), ( pT]) into (|17D and differentiating the equation ob- 
tained with respect to \j we get Aqx-x^ = 0,i,j = 1,2, hence 

Ao = XiRiit) + A2i?2(t) + Roit), (28) 

where i?i, i?2, -Ro are arbitrary smooth functions. 

With these results we can split equations ([TT|) and (|12|) by Ai,A2 thus 
obtaining a system of four nonlinear PDEs for the three functions uii,ui2,Q 

Q0J2t + yQ0J2x = iyy + kx)QuJ2y + 2vQyUJ2y + iyQuj2yy, (29) 

Qi + (5i?o + QLoiujit + yuji^) + yQx = vQ + {yy + kx)Qy 

+^Qyy + jyQSoujly, (30) 

Ri + Li{uju + yi^ix) = ^Siujjy, (31) 

R2 + L2{uju + yui^) = yS2Ujly. (32) 

Making an equivalence transformation (P) with appropriately chosen func- 
tions we can put Lq = 0, i?o = 0. Next, due to the requirement (|^) S1S2 7^ 0. 
There are two inequivalent cases L2 = and L2 7^ 0. Since they are 
handled in a similar way, we consider in detail the case L2 = only. In view 

13 



of (^ Li does not vanish. Choosing appropriately the functions /i, /2 in (H) 
we can put Li = 1,82 = ±1 in formulae (p9D -(p2|). Integrating (|32D with 
account of ( PT| ) yields for U2 

uj2 = R{t)y + F{t,x), R{t)^0, (33) 

where R, F are arbitrary smooth functions and R2 = ±z/i?^. 

Differentiating ([31| ) twice with respect to ?/ and taking into account (^) 
we arrive at the equation Siuj^uj^ = 0, therefore 

oi = CiCt;2 + C2, 

where Ci 7^ 0, C2 are arbitrary constants. Next, integrating (pi]) we obtain 
for cji, F(t,x) 

u;, = uC,{R'x + P{t))-JR^{t)dt, 

F{t, x)=3R + R-^P{t) - C{^C2, 

where P{t) is an arbitrary smooth function. 

Hence we conclude that the corresponding solution with separated vari- 
ables reads as 

u = Q{t, X, y) exp{Ai / i?i(t)dt + A2 / -R2 Wdt} exp{AiCo'i}(y92(co'2) 
= Q(t,x,y)exp{A2 /i?2(t)dt}exp{Ai(z/Ci(i?='x + P(t)))}^2(cu2). 



Thus, the function Ri{t) does not enter into the solution with separated 
variables and, therefore, we can put i?i = in (|3^). Furthermore, within an 



equivalence transformation (^ we can choose Ci = v ^, C2 = 0, thus getting 

oji = R{tfx + P{t), (35) 

U2 = Rit)y + 3R{t)x + PR{t)-^. (36) 

Provided L2 7^ 0, the forms of the functions u!i,u!2 is the same as those 
given in (©,(©. 

Inserting (^) and (^6]) into (pQ]) and integrating by y we get the form of 
the factor Q(t, x, y) 



\(AR-vR\ 2 3R-kR y d ( P^ 

^ = '"'''[[-^iJR-jy ^ ^l^jR-'^y ^ 2iMJt[-R\ 

+M(t,x)l. 



14 



(37) 



Substituting ( pTj ) into pO| ) we come to the following relation 

where we use the notation 

Differentiating ( pSD three times with respect to y yields Sq^^^^^,^^ = 0, 
therefore 

^0 = C1L1J2 + C2UJ2 + C3, 

where Ci, (72,6*3 are arbitrary constants. Next, differentiating ( pSD with re- 
spect to y twice and with respect to x once we get M^xx = 0, or 

M = Mi{t)x^ + M2{t)x + M3(t), 

where Mi, M2, M3 are arbitrary smooth functions. 

Finally, inserting the obtained expressions for 5*0, M into (|38|) and split- 
ting by the variables x, y we come to the following system of ODEs: 

5 = 27p+-C,fl'-^+5. (39) 

M2 = -— Z + —Z + z/Cai?^ + 2uRPCi, (42) 

o p 

M2 = — -Z + 3iyR^RC2 + QuRPCu (43) 

2h'R 

^3 = ^ + 2^ + -^Z^ + ^^71^ + UC2P + lyCsR^. (44) 

15 



Differentiating (|iO| ) with respect to t and subtracting the resulting equa- 
tion from (41) yields the fourth-order ODE for the function R 



R^^^^ ^RR ^R^ ^^R^R , 2^ f^^s k 



3 



6^ + 2^ - 10^-^ + Au'^CiRR^ + ^ = 0. 



R R^ R^ i?3 ' 3 

Reducing the order of the above ODE with the help of equation (|39D and its 
first- and second-order differential consequences we arrive at the following 
relation: 

47/2 4 U2 2u U2 

—C^rS = — + — - — -— f45^ 

25 ^ 100 25 25 9 ■ ^ ^ 

If in (H) Ci ^ 0, then in view of (|3|) A; = 0. Provided, Ci = 0, fc is a 
root of the quadratic equation 

64fc2 + 36u^k - 9z/2 = 0, 

hence k = Su^/IG or k = — 3z/2/4. 

Thus system of ODEs (p9D-(^) is consistent only if the parameter k takes 
one of three values 0, 3i/^/16, — 3z/^/4. Consequently, KE (^ has solutions 
with separated variables in the case considered (i.e., provided the system (^) 
takes the form (^) only for the values of the parameter k given previously. 
This provides the proof of the first part of Theorem 1. 

We examine the three possible cases 0, 3z/2/16, — 3z/2/4 separately. 

Case 1. For k = 0. Then the equality R{t) = ±2-^/^S^^^^ = const holds. We 
denote this constant as r. Next, it follows from (|43|) that M2 = m = constant. 
In view of these facts we get from ( ^21) ODE for P(t) 



- P +u^P + 2ur^{uS2r^ - m) = 

which general solutions reads 

P{t) = C^e'"' + Cse-'^* + 2r%mu-^ - S2r^)t + Cq, (46) 

where 04^,0^,0^ are arbitrary constants. 

A direct check shows that applying finite transformations from the sym- 
metry group admitted by KE under fc = to the obtained solution with 
separated variables (|^), (p5|), (p6D, (|i6| ) we can cancel P(t). 
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Scaling when necessary u!i,u!2 in (p5D, ( p6| ) we can choose r = 1. Hence 
we get the equahty Ci = |. Summing up we conclude that the following 
relations hold: 

Q = exp I -— + z/Cax + z/ fCg + -j tj , ui = x, uj2 = y, 



<po = u (Xi - f C3 + 2 



^0, 



(pi = V{\2 - 6*2)931, <^2 = I ^ + A2a;2 + Ai - - I (/32- 

Then the corresponding solution with separated variables is 

u = (f2 exp -— + u{Xit + X2X) 1 . 

Consequently, the constants C2 and C3 + | do not enter the final form of the 
solution with separated variables. This means that we can put C2 = and 

C3 = -2- 

Thus we have proved the validity of the first part of Theorem 5. 
Cases 2,3. For A; = ^ or A; = — ^. In these cases we get from ( pUD 

R (rV 2 

r-\r) =-"• 

where 

4 lb 

2' 4 

Integrating the above ODEs yields 

R{t) = (Ci sinh at + C2 cosh at)^^, 

where Ci,C2 are arbitrary constants. 

Using shifts with respect to t and the equivalence transformation (||) we 
get the four inequivalent forms of the function R{t) 

Rii) = — , R(t) = — , R(t) = expjiat}. 

^ ^ coshat' ^ ^ sinhat' ^ ^ t'l j 
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Comparing (^) and the first-order differential consequence of ( ^31) yields 
the second-order ODE for Z{t) = R-\d/dt){P/R^) 

= 0. (47) 

The general solution of this equation has the following structure: 

Z{t) = CiZi{t)+C2Z2{t), 

where Ci, C2 are integration constants. Hence, we conclude that the function 
P{t) is of the form 

P{t) = C,P,{t) + C2P2it) + CsP^it) + C,P,{t), (48) 

where C3, C4 are integration constants. 

On the other hand, if we apply to the solution with separated variables 
(^, (|35|), (p6D with P{t) = finite transformations from the symmetry group 
of KE under k = ^ or k = —^, then we get an equivalent solution with 
separated variables such that P{t) is of the form 

P{t) = C[P[{t) + C',P^{t) + C'.P^it) + C'.Piit). (49) 

Here C(, . . . , C4 are arbitrary constants and the functions P[{t), . . . ,P^{t) 
are linearly independent. Hence, we conclude that due to the theorem on 
existence and uniqueness of the Cauchy problem for a fourth-order ODE 
(^71) (considered as an equation for the function P{t)) the expressions on 
the right-hand sides of (^) and ( ^Uj ) coincides within the choice of constants 
Ci,C[,i = 1,...,4. Consequently, without loss of generality we can put 
P{t) = in formulae (^), (^). 

Using the reasonings analogous to those of Case 1 we can put r = 1, C2 = 
0, C3 = 0. The second part of Theorem 6 is proved. 

A similar analysis of the separability of KE into three ODEs (|^ yields 
the proofs of the remaining assertions from Section 2. 



4 Exact solutions 

Remarkably, for the equation under study it is possible to give a complete 
account of solutions with separated variables. For the case when KE sepa- 
rates into three first-order ODEs (|^, we get the following family of its exact 



solutions: 

/ /■ //iAi + /2A2y ,, , , fiy-fix f2y-f2X 
u = exp< u / -. : — at + Ai^ : h A2-- ■■ \- 

\ J \ /2/1 - /1/2 ; /2/1 - /1/2 /2/1 - /1/2 



4^/2/1-/1/2 4;^ 2z/ V/2/1-/1/2 

+ I -^^^^i^M - ^) x^ - llnlAA - A/2I + -tl 
\4z. /2/1 - /1/2 4J 2 '^^'^ '^''' 2/' 

where /c, /i(t), /2(^) are given by the corresponding formulae from Theorems 
2-5. 

Next, for the case when KE separates into three ODEs of the form (|^) 
we obtain the following families of its exact solutions: 

1. A; = (this case has been considered in Theorem 5) 

M = exp ( -— + u{\it + Aax) j Dy2^^^ (y + 2A2), 

where D^ is the parabolic cylinder function. 

2. k = 3z/Vl6 or A; = -Su^/i 



u = exp<j z/Ai / RMt + v\2R^x ( :^ - 7 ) Z/^ + ^ ( 3— - A; j x?/ 



+ (-11 + ^-1)^^ n' + ^'"«|w*^ + «^) 



2 
■^ V3T2' 



+Ai}5Zi -— (A2(% + 3i?x) + Ai)2 



where i? is given by the corresponding formula from Theorem 6 and Zi is 
the cylindric function. 

Note that the above obtained families of exact solutions of KE contain 
two continuous parameters Ai, A2. These parameters have the meaning of 
eigenvalues of two commuting symmetry operators of KE, while the cor- 
responding solution with separated variables is the eigenfunction of these 
operators. Provided some appropriate boundary and initial conditions are 
imposed, the parameters become discrete and thus we get a basis for expand- 
ing sufficiently smooth solutions of KE into series. 
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5 Conclusions 

It is a remarkable feature of the Kramers equation (|^) that a classical problem 
of variable separation can be solved in full generality. The results obtained on 
this way are in good correspondence with the ones on symmetry classification 
of KEs of the form @. As follows from the papers |^, ^, the cases k = 
Su'^/IG, k = — 3z/^/4 are distinguished by the fact that the corresponding 
KEs (0) admit the most extensive symmetry groups. For these choices of 
k KE (j^) is invariant with respect to eight-parameter Lie transformation 
groups, while for all other values of the parameter k the maximal group is 
six-parameter. 
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